Dirac electrons in a Kronig-Penney potential: 
dispersion relation and transmission periodic in the strength of the barriers 
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The transmission T and conductance G through one or multiple one-dimensional, (^-function 
barriers of two-dimensional fermions with a linear energy spectrum are studied. T and G are 
periodic functions of the strength P of the 5-function barrier V(x,y)/hvF = P8(x). The dispersion 
relation of a Kronig-Penney (KP) model of a superlattice is also a periodic function of P and causes 
collimation of an incident electron beam for P — 2nn and n integer. For a KP superlattice with 
alternating sign of the height of the barriers the Dirac point becomes a Dirac line for P — (n+l/2)n. 
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The study of particle motion in periodic potentials is 
at the heart of condensed matter physics and it is usually 
assumed that the energy spectrum is parabolic. One of 
the earliest examples is the well-known, one-dimensional 
(ID) Kronig-Penney (KP) modelP that consists of an 
infinite succession of very thin (W — * 0) and very high 
(Vq — > oo) barriers, referred to as (^-function barriers, 
but such that their product P cx WVq remains constant. 
This results in minibands in the electron spectrum. 

One may wonder though how such results are mod- 
ified if the energy is linear in wave vector. Such a 
spectrum occurs for relativistic electrons with energy 
E = hep >> Eq — moc 2 , where c is the speed of light and 
mo the bare electron rest mass. Even without neglecting 
Eq a strict ID Dirac KP model was considered for rela- 
tivistic quarks^. It is also known that electrons can trans- 
mit perfectly, upon normal incidence, through arbitrar- 
ily wide and high barriers, referred to as Klein paradox 
or Klein tunneling^. With the discovery of graphene 4 , 
a one-atom thick layer of carbon atoms, another system 
became available in which particles (electrons) moving in 
two dimensions, have a linear spectrum, E = hvpk, with 
k = (k x ,k y ) the wave vector. Importantly, carriers in 
graphene behave as chiral, massless fermions described 
by Dirac 's equation without the mass term, and move 
with the Fermi velocity vf ~ c/300. There is a wealth of 
exceptional properties of graphene, see e.g., Refs 5 . 

Because the carriers in graphene move in two di- 
mensions, tunneling through barriers is inherently two- 
dimensional (2D) and depends on the direction of the 
incident electron beam even in the absence of a mag- 
netic field. Many authors, including ourselves, have 
studied this tunneling, through single, multiple barriers, 
and superlattices^l Surprisingly, tunneling through 8- 
function barriers has received very little attention^ and 
we are not aware of any Dirac KP model for a superlattice 
in graphene. An interesting development was the appli- 
cation of periodic potentials to graphene that turned it 
into a self-collimating material despite the rather unusu- 
ally high potentials usecP. 

Motivated by all these results and the absence of a 
systematic treatment of KP barriers or superlattices, we 



study in this work the transmission through such struc- 
tures as well as the dispersion relation of a KP super- 
lattice. Although the model may appear a bit unrealis- 
tic, since a relatively smooth potential is needed to de- 
scribe the carrier dynamics by the Dirac equation^, its 
simplicity is attractive and elucidates certain symmetry 
properties of the spectrum. Furthermore one can real- 
ize the model by using a potential which is smooth on 
the scale of the atomic distance while remaining imidiate 
compared to the typical electron wavelength. The un- 
expected results mentioned in the abstract are in sharp 
contrast with those for carriers with a parabolic energy 
spectrum described by the Schrodingcr equation. We will 
use graphene as an example but the results apply to any 
2D system with a linear-in-wave-vector spectrum and a 
two-component spinor. 



I. TRANSMISSION THROUGH A 5-FUNCTION 
BARRIER 

We describe the electronic structure of an infinitely 
large flat graphene flake in single valley approximation by 
the zero-mass Dirac equation and consider solutions with 
energy and wave vector near the K point. The Hamilto- 
nian is Ti. — vpa ■ p+Vf^V with p the momentum operator 
and ¥ the 2x2 unit matrix. In the presence of a ID po- 
tential V(x) the equation (H — E)ip — admits solutions 
of the form i\i{x)e lkyV where 

4>W = LV) ^)=(_ se W)e- iAx , (1) 

with tan0 = k y /X, s = sign{e — u(x)) , A = [(e — u{x)) 2 — 
ky] 1 / 2 , e = E/vfH, and u(x) = V(x) /vph; e and u(x) are 
in units of inverse length. As usual, we approximate a 8- 
function barrier with a very thin and very high barrier, of 
width W(-^ 0) and height Vq(—> oo), but keep constant 
the dimensionlcss product P — WVo/Hvp which we call 
its strength. Referring to Eq. (fll) and Fig. [l] the wave 
function in each of the regions (l)-(3) can be written as 
a superposition of the eigenstates of Eq. ([!} tp n (x) = 
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Figure 1: Left: A ID potential barrier of height Vo and width 
W. Right: wave vector k = [k x , k y ) geometry for an electron, 
with energy < E < Vo, impinging on a square potential 
barrier (grey area). 

Gn-M- n (x)A n ,n = 1,2,3, with coefficients A — (A,B) T 
and 

Q=(J# _ se %), M(,)=( e 7 e JL). (2) 

Matching the wave function at the interfaces x — and 
x = W gives the transfer matrix JV = Q~ X SQ in terms 
of the matrix S which relates the wave function in front 
of the barrier to the one after it in the manner ipi (0) = 
Sy>3(0+). The result is 

We notice that S is a periodic function of P and that 
S = ±F for P = mr. This is a special situation in which 
the two pseudo-spin components of the wave function 
do not mix. Later on we will see that this periodicity 
appears in the transmission through a <5-function barrier 
and in the dispersion relation of the KP model. 

With the elements of Af denoted by n,j , the transmis- 
sion is T = \t\ 2 = l/|nn| 2 . The explicit result is 

T= l/[l + sin 2 Ptan 2 <?(>] , (4) 

and coincides with the formula for transmission, found 
in RefP, in the limit of 5-function barriers. Obviously, 
T and R = 1 — T are periodic functions of P, that is, 
X(P + nn,<f>) = X(P,d>) for n integer and X = T,R. 
In addition, from Eq. Q we find that T(P, 0) has the 
following properties: 

1) T(P, 0) = T(tt -P,4>)= T(tt + P, 0) , 

2) T(nn,ct>) = l, 

3) T(tt/2,0) = cos 2 0, 

4) T(P, 4> = 0) = 1, T « 1 for <t> w k v < k x , 

5) T(P, ±7r/2) = 0, T « for cf> w ±tt/2 ^ k y ■> k x (5) 

These results are very different from those of the non- 
relativistic case where T is a decreasing function of P. 
A contour plot of the transmission is shown in Fig. 
2(a). This figure shows clearly the symmetry properties 
T(P, 0) = T(P, -0), and T(n -P,</>) = T{P, 0). 



A. Conductance 

The two-terminal conductance is G = 
G J^ 2 /2 T(P,(f>) cos (f>d<l), with G = 2E F L y e 2 /(v F h 2 ) 
and L y the width of the system. Using Eq. (4) for 
T(P, <j)) the resulting G is periodic in P and given by 

G/G = 2[1 - artanh(cosP)sinPtanP]/cos 2 P. (6) 

For one period, G is shown in Fig. 2(b); its minimum 
value is 4/3 and its maximum one 2. 
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Figure 2: Left: Transmission T through a 5-function barrier 
vs its strength P and the angle of incidence (j> (tan <j> — k y /k x ). 
Only one period is shown. Right: Conductance G/Go vs P. 



II. TRANSMISSION THROUGH TWO 
(5-FUNCTION BARRIERS 

We consider two barriers separated by a distance L 
characterized by the potential V(x,y)/hvp = Pi5(x) + 
p2<5(a; — L), with strengths P\ 2 and introduce the di- 
mensionless variables e — > eL, k y — > k y L, uq — > woi, and 
x — > x/L. Due to space limitations we treat only the 
cases of parallel and antiparallel (5-function barriers with 
the same strength | J\| = |P 2 |. 

Parallel 5 -function barriers. Thi s is a model system for 
a resonant tunneling structure^' and also for a Fabry- 
Perot interferometer whose resonances were recently in- 
vestigated experimentalljffS. The transmission is given 

by 

T = [l+tan 2 0(cos k x sin 2P-2s sin k x sin 2 Pj cos 0) 2 ] 1 , 

(7) 

with s = sign(e). The properties of T{P,(j),k x ) are 
identical to those for a single barrier except for prop- 
erty 3) and property 1) which must be replaced by 
T(P,<j),k x ) = T(P + rm,(f>,k x ). In Fig. 3(a) we show 
T(P, (j), k x ) through two barriers with P = 7r/10. 

Antiparallel 5 -function barriers. We now consider two 
parallel (^-functions but with opposite sign, P = P x = 
— Pi. The result for the transmission is 

T= [cos 2 /c 2 . + sin 2 fc a; (l-sin 2 (/)cos2P) 2 /cos 4 0] _1 . (8) 

In Fig. 3(b) we show the transmission through two op- 
posite barriers for P = 7r/10. The symmetry properties 
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Figure 3: (a) Transmission through two parallel 5-function 
barriers, as a function of the wave vector component k x and 
angle of incidence <f>, for P — it/ 10. The yellow solid curve 
shows the contour with a transmission value of 0.75 for P = 
7r/2. (b) As in (a) for two antiparallel 5-function barriers. 



of T(P, cf>, k x ) for a single barrier again hold, except for 
the value of T(w/2,(f),k x ), see property 3). In addition, 
we now have T(P, <j), k x ) = T(P, </>, -k x ). 

The periodicity in the transmission is also present in 
the conductance G. We show G in Fig. 4(a) for parallel 
and in Fig. 4(b) for antiparallel (^-function barriers. 
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leads to (A = [e 2 - fc 2 ] 1 / 2 ) 

cos k x = cos P cos A + (e/ A) sin A sin P . (9) 

The solution of Eq. ([9| gives the dispersion which is 
periodic in P and the spectrum is shown in Fig. 5 for 
P = it/2. Further, Eq. ^ is mapped, for k y = 0, 
directly onto that for strictly ID fermions 2 and gives the 
spectrum 



e = P ± k x + 2nir , 



(10) 



with n an integer. 

Equation (9) contrasts very sharply with that for 2D 
electrons with a parabolic spectrum in a ID KP potential 
which, with A' = [2 fie — k 2 ] 1 / 2 and fi = mvpL/h, reads 

cos k x = cos A' + (fiP/X') sin A' ; (11) 
the resulting dispersion relation is not periodic in P. 
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Figure 4: (a) Conductance G(Go) vs e for several strengths 
P through two parallel 5-function barriers, (b) As in (a) for 
two antiparallel barriers. 



Figure 5: (a): Energy bands, near the Fermi level, close to 
the K point in the KP model, for P — tt/2. The valence 
band touches the conduction band at e = P = tt/2. For large 
k y the valence band becomes flat, (b) and (c): projections, 
respectively, of the valence and conduction bands onto the 
(k x ,k y ) plane. 



III. KRONIG-PENNEY MODEL 



For an infinite number of periodic square barriers, one 
can tune the system into a self-collimating material 9 . For 
special values of V , W, and L it was found in RefP 
that the dispersion relation near the K point is almost 
linear in k x and constant along k y . The system thus 
behaves as a ID metal. We look for similar results using 
(^-function barriers. Within the KP model we replace the 
square barriers by (^-function barriers, characterized by 
V(x,y)/hvF = X)jL_oo PS(x — jL). The resulting wave 
function is a Bloch function and the transfer matrix Af 
pertinent to these barriers leads to ^(l) = e lkx ij){/S) and 
A\ = AfA.2, with k x the Bloch wave vector. From these 
conditions we can extract the relation e~' lkl M(\)A2 = 
NA 2 , with M{x) given by Eq. (pi). Then setting the 
determinant of the coefficients in Ai — (A B) T equal to 
zero and using the transfer matrix for a <5-function barrier 



A. Properties of the spectrum 

Since the dispersion relation is periodic in P, with pe- 
riod 2tt, it is sufficient to study it only for < P < 2tt. 
For particular values of P we find: 

1) P = 2tt -> e = ±[fc 2 + (k x + 2nir) 2 ] 1 ' 2 , 

2) P = it - e = ±[fc 2 + (k x + (2n + l)^) 2 ] 1/2 , 

3) P = 7r/2(37r/2) -> cosfc,, = +(-)esinA/A. (12) 

In limiting cases we are able to obtain explicit expres- 
sions for E — E(k x ,k y ). We expand the dispersion re- 
lation for small k y and e — P. The resulting quadratic 
equation for e is solved by 

e± w P ± [4sm 2 (k x /2) + (k 2 y /P 2 ) sin 2 P] 1/2 . (13) 
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For small k x we can replace the term 4sin 2 (fc :r /2) by k 2 . 
Notice that for k x — we find e± w P±k y sin P/P which 
is a linear spectrum with a reduced velocity. For fc y = 0, 
we have e± — P ± 2sin(|fea;|/2), which is linear for small 
fe x but posseses a typical band shape for large k x w 7r. 
For small fe^ > fcj, sin P/ P we have 



e± «P± 2fc 2 



fc^ sin' 



P/P ; 



(14) 



For P 3> 1, £ is highly anisotropic and nearly flat vs k y . 

Relation to the spectrum of a square superlattice. We 
now look whether we can find an energy spectrum sim- 
ilar to that of Ref.S pertinent to square barriers, with 
height Vq = 720 meV, and width w = 5 nm, and 
unit-cell length L = 10 nm. In our units these val- 
ues correspond to P = VoW/vph — 2tt and lead to: 
e = ±[ky + (k x + 2m:) 2 ] 1 / 2 . Since the Fermi level in 
these units is ep = 2ir — P, we look for the spectrum 
near the value e = [k 2 + (±\k x \ + 2-k) 2 ] 1 / 2 . Although 
these bands seem to fullfil our demands because the dis- 
persion looks rather flat in the k y direction, the concern 
is that we would obtain the same dispersion for P — ► 
and £f —* 2-7T. But this can be obtained by folding the 
cone-like dispersion of graphene and results simply from 
working in the reduced-zone scheme. Consequently no 
new fundamental physics should be attached to it. Fur- 
ther, from this correspondence we expect and found that 
for square barriers with P = 27m, the situation is more 
favorable for the occurrence of collimation. It follows 
that the collimation effect is also obtainable for barriers 
that are lower than the unusually high ones of Ref. 9 if 
one uses longer unit-cell periods. 



IV. EXTENDED KRONIG-PENNEY MODEL 

The square barriers are replaced by alternating-in-sign 
5-function barriers. The unit cell of the periodic potential 
contains one such barrier up, at x = 0, followed by a 
barrier down, at x — 1/2. The resulting transfer matrix 
leads to 



cos k T 



cos A - {2k 2 /A 2 ) sin 2 (A/2) sin 2 P , 



(15) 



where tan</> = k y /X. From Eq. (15) we deduce that the 
dispersion is periodic in P, with period 7r, and has the 
following properties: 

1) it is invariant for e — > —e and P — > n — P , 

2) P = nir -> £ = ±[fc 2 + (k x + 2mr) 2 ] 1/2 , 

3) P = tt/2^ (e,k x ,ky) = (0,0, fey)- (16) 

In Fig. 6(a) we show the spectrum for P = ir/2. As 
seen, it is almost independent of k y for small energies 
while in the k x direction the bands are linear: this is an 
advantageous situation for self-collimation. For k y = 
we obtain the linear spectrum 




Figure 6: (Color online) (a) Conduction and valence bands, 
near the K point, in the extended KP model, for P — n/2. 
The bands form a "cross" in the (E, k x ) plane and the Dirac 
point has become a Dirac line. The dispersion is nearly inde- 
pendent of k y . (b) As in (a) for P — n/4. 



with the Dirac point at e = 
explicit expression for k x « 0. 



0. We can also find an 



Solving Eq. (15 1 gives 



±\k y cos P| 



(18) 



Then the group velocity v y oc de/dk y becomes small if 
P w 7r/2 + nir. Fig. 6(b) shows the energy spectrum 
for P = 7r/4, the Dirac cone becomes anisotropic, as the 
spectrum flattens in the k y direction. 

We now consider the case where k x and k y are nonzero. 



If £ <C 1 then the rhs of ( 15 ) can be expanded in e. This 



leads to a quadratic equation for e with solutions 



±\ky 

P=tt/2 



cosh k y — cos k x — f 



1/2 



(fey/2) cos 2 Psinh k y + f 
±\k y \sm(\k x \/2)/ smh(\k y \/2) , (19) 



where / = 2 sin Psinh (k y /2). For k y = we find the 
result given by Eq. (131, that is, e± = ±2sin(|fe x |/2) 



±\k. 



2nn . 



(17) 



which is linear for small k x 



V. CONCLUSIONS 

In summary, we studied the transmission and conduc- 
tance of fermions, with energy linear in wave vector, 
through one and two (5-function barriers and the energy 
spectrum of a KP superlattice. For very high (Vq — > oo) 
and very thin (W — ► 0) barriers we showed that they 
are periodic functions of their strength P = WVo/Kvp 
where vf is the Fermi velocity. Further, we showed that 
a KP superlattice has an energy spectrum that is a pe- 
riodic function of P, which is in sharp contrast with 
that obtained from the Schrodinger equation. An im- 
portant consequence of that is collimation of an incident 
electron beairP that here occurs for P = 27m with n 
an integer. We also obtained various explicit but ap- 
proximate dispersion relations, e.g., for small wave vec- 
tor k — (k Xl k y ). Given the intense research activity in 
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graphene and the very recent experimental verification Policy (IAP), and the Canadian NSERC Grant No. 
of Klein tunneling^, we expect that this periodic depen- OGP0121756. 
dence on the strength P will be tested in a near future. 
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